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^ ■ Abstract 
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ON ■ We present simple calculations which show that the incommensurability upon 



doping and the width of the magnetically ordered phase in Mott-Hubbard 



a 

i-rt ■ insulators depend strongly on the location of the hole/electron pockets in the 

a: 

Brillouin zone. For LaSrCuO systems, we found the pockets at (±7r/2, ±7r/2), 



in which case the corrections to the antiferromagnetic spin stiffness rapidly 
grow with doping and destroy commensurate spin ordering already at a very 
small doping. On the other hand, in NdCeCuO, the hole pockets are located 
at (tt, 0) and the symmetry related points, in which case the corrections to 
the stiffness scale linearly with the density of carriers and do not destroy 
commensurate spin ordering. For YBCO systems, the situation is less certain, 
but our results favor hole pockets at (tt/2,tt/2). We also briefly discuss the 
tendency towards phase separation. 
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I. INTRODUCTION 

The intense interest in understanding the properties of high temperature superconductors 
initiated theoretical research on the behavior of antiferromagnetic insulators upon doping. 
The parent compounds of high-T c materials are well described as Heisenberg antiferromag- 
nets. Upon hole doping, long-range antiferromagnetism rapidly disappears and the systems 
eventually become metallic superconductors. The same transformation occurs in electron- 
doped materials but at substantially larger doping concentrations. The behavior of the 
antiferromagnetic insulator upon doping attracted a lot of interest over the past few years 
as the exchange of antiferromagnetic paramagnons is at least one of the relevant pairing 
interactions between holes Q. There are several fundamental issues related to doped anti- 
ferromagnets, one of which is whether magnetic correlations remain peaked at Qq = (n, n) 
upon doping, or shift to incommensurate momenta. Shraiman and Siggia first pointed out 
that if the dispersion of vacancies has a minimum at (±7r/2, ±7r/2), then the dopants intro- 
duced into commensurate nearest-neighbor antiferromagnet give rise to a long-range dipolar 
distortion of the staggered magnetization which may lead to a spiral spin configuration. 

It has been recently argued |3|-[5| that the values of hopping integrals in electron-and 
hole doped 214 materials are nearly the same; t ~ OAev, t' ~ — 0.2t, where t' is the nearest- 
neighbor hopping. Despite this, the incommensurability upon doping has been found only 
in LaSrCuO compounds || while the dynamical structure factor in the doped NdCeCuO 
remains peaked at (it, it) f|0]. Moreover, experimentally, long-range magnetic order in 
LaSrCuO, disappears already at 3 — 4% doping, while in NdCeCuO it survives up to 12% 
doping. 

The goal of the present paper is to show that the contrasting magnetic dynamics in 
the two 214 systems near half filling is related to a different location of the hole pockets. 
Namely, we will argue that in La— based materials, the pockets are located at (±7r/2, ±7r/2) 
while in Nd— based materials doped electrons occupy pockets centered around (0, ±7r) and 
(±7r, 0). This, as we show below, gives rise to a completely different spin dynamics in the 



two materials: the commensurate (tt, tt) state rapidly becomes unstable in LaSrCuO, but 
survives in NdCeCuO. 

We will also discuss the location of the hole pockets in YBa2Cu^O^ +x where the next- 
nearest neighbor hopping amplitude is relatively large, t' ~ — 0.5£. Our results for the 
Hubbard model, show that for t'/t ~ —0.5 and J/t ~ 0.4, the mean-field hole dispersion 
is nearly degenerate along k x = k y and has a flat minimum at (7T,7r). However, self-energy 
corrections still favor pockets at (tt/2, tt/2) and are likely to overshadow the small difference 
between the mean- field quasiparticle energies at (tt/2, tt/2) and (tt, tt). In this situation, 
the spin dynamics of LaSrCuO and YBaCuO near half-filling are nearly identical, and 
differ only in the metallic phase where the Fermi surface is large, and in case of YBaCuO, 
is centered at (tt,tt). Notice, however, that cluster calculations for the t — t' — J model 
reported hole pockets at (tt, tt) for the same ratios of parameters P,|lT|. If it is actually the 



case for YBaCuO, the spin dynamics right near half-filling will be very similar to that in 
the electron-doped materials (see below). 

The bulk of our consideration is presented in the next section. We will first briefly 
review the spin-density wave theory for the Hubbard model near half-filling, then will find 
the location of the hole pockets for the hole- and electron-doped materials, and next show 
how the different location of the hole pockets give rise to a contrasting magnetic behavior 
near half-filling. Finally, we discuss the tendency towards domain wall formation upon 
doping. Our conclusions are presented in Sec.p|. 

II. HUBBARD MODEL WITH NEXT-NEAREST-NEIGHBOR HOPPING 

We consider the one-band Hubbard model given by 

Here j and j' label the nearest and the next-nearest neighbors, respectively, and n = c^c is 
the particle density. We will use the spin-density-wave formalism [12[] which, as has been 



shown in a number of papers fll2|-(18|1 , is a good starting point for the calculations close to 



half-filling. Below we mostly restrict with the simplest mean-field calculations. This last 
restriction can be formally justified if one extends the Hubbard model to a large number of 
orbitals at a given site, n c = 2S, [[HJ and restrict with the leading term in 1/5* expansion. 
This mean-field theory is meaningless for the nearest-neighbor Hubbard model because of 
the accidental degeneracy in the hole spectrum which is lifted only by 1/S corrections (see 
the discussion below). However, the nonzero t' eliminates the accidental degeneracy already 
at the mean-field level. In this situation, we expect that the corrections to the mean-field 
results renormalize the parameters of the model, which will be important for our analysis of 
YBaCuO, but do not give rise to any new physics of the insulating phase. 

A. SDW theory at half-filling 

We now briefly discuss the key points of the SDW formalism at half-filling. This formal- 
ism has been applied several times to the t' — model. First, we assume that at half-filling, 
the 2D Hubbard model has a commensurate antiferromagnetic ground state. This im- 
plies that, e.g., z component of the spin-density operator, S(q) = (l/2)J2k a k+q,a&a,i3 a k,i3, 
has a nonzero expectation value at q = Qq. We then use the relation (J2k a k+Q i a k,i) = 
~~ (J2k a 'k+Q, i a k,i) = (Sz)i to decouple the quartic term in (|l|). After decoupling, the quadratic 
Hamiltonian takes the following form 

'HMT = Yl e k( a ka a ka + a k+Q Q a a k+Q a) + 
k 

/ , 6 k ( a kcr a kcr — a k+Q cr a k+Qocr) ~ 

k 

Y!^sgn(a) (a ] k(T a k+QQ(T + a\ +Qo(7 a k(7 ) (2) 

k 

Primes to the summation signs indicate that the summation is over the reduced Brillouin 

zone. We introduced A = U(S Z ), e k = —2t(cosk x + cosk y ) — At' cosk x cosk y , e + = (e k + 

e k+ Q )/2 = —4t'cosk x cosk y , e~ = (e k — e k+ Q )/2 = —2t(cosk x + cosk y ). The next step is 

the diagonalization of the quadratic form by a Bogolyubov transformation 



0>k,a — UkCka + Vkdka 

a>k+Q ,cr = sgn(a)(u k d k a - v k c ka ) 



(3) 



Applying this transformation to @, we observe that the first term with the density of 
quasiparticles, indeed, does not depend on u k and v k , because the transformation conserves 
the total density. The Bogolyubov coefficients then appear only in the last two terms which 
do not depend on t' . As a result, the expressions for u k and v k remain the same as in the 
£' = model IB3I: 
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where El 



A 2 + fer) 2 . 



After the diagonalization, Eq.(^) takes the form 

'Hmt = Yl E l c ku c ku - E^d ka d ka 

k 

where 



(5) 



E c = EZ + e+, E° 



El 



(6) 



For U ^> t, which is implicit in our approach, we can expand under the square root and 
obtain E c ' = A + J(cosk x + cosk y ) 2 =p At' cos k x cos k y , where J = At 2 /U. We will refer 
to the quasiparticles described by c and d operators as conduction and valence fermions, 
respectively. At half-filling, valence states are occupied and conduction states are empty. 
Accordingly, the self-consistency condition on (S z ) takes a simple form 



TT ~ ^ TP- 



u 



/, E k 



(7) 



At large U, we obtain, as usual, A = U/2, or (S z ) ~ 1/2. 



B. Finite density of holes 

We now discuss what happens at small but finite doping when the chemical potential 
moves into the valence band. First, we discuss the shape of the hole Fermi surface. As we said 



above, at t' = 0, the dispersion of valence fermions, Eq.@, is degenerate along the boundary 
of the magnetic Brillouin zone (k x ± k y — ±vr), where E d = A. This degeneracy, however, 



is not related to any kind of symmetry and is removed by self-energy corrections fT5| , |TT|| , 
with the result that the actual band minima in the nearest-neighbor Hubbard model are at 
(±7r/2, ±7r/2). This agrees with the numerical |2"4| , |2T)| , |2"o| and variational [2B] studies of the 
Hubbard and t — J models. 

Our first observation for t — t' — U model is that at finite t', the degeneracy is removed 
at the mean-field level. Indeed, a simple inspection of Eq.@ shows that the mean-field 
dispersion has a minimum at (±7r/2, ±7r/2) if t' is negative and smaller than J . For \t'\ > J, 
the minimum of the hole dispersion is at k = (it, it) (or (0, 0)). Finally, if t' is positive, which 
is probably not the case for cuprates, the minimum of E d is at (n, 0) and symmetry related 
points |20|| . Now, for both 214 compounds, t' is negative and relatively small: \t'\ ~ 0.07eu, 
which is smaller than the exchange integral J ~ 0.13ei> [pl| -p3[]. Accordingly, we expect 
that upon doping, holes in La— based compounds form pockets around (tt/2, n/2) and the 
symmetry related points. Near (tc/2,tt/2), one can expand E d and obtain 

k 2 t kl 

where m± = 1/4(J — |t'|), mn = l/4|t'|. Notice that numerically, the effects due to £', 
even for \t' /t\ ~ 0.2 are likely to dominate over the effects due to self-energy corrections in 
nearest-neighbor model. Thus, for t/J = 2, the difference AE = E d (7c, 0) — E d (ir/2, vr/2) is 
AE = 4\t'\ ~ 0.8t due to t' and about 0.25t due to quantum fluctuations as was obtained 



in 1/S expansion for the Hubbard model JT7j and in numerical |25| and variational |26 



calculations for the t — J model. The same is also true for the inverse effective mass, 1/mf 
for the same ratio t/J, the contribution to 1/toii due to t' is about four times larger than 
due to quantum fluctuations. Also notice that for the parameters chosen for LaSrCuO, the 
two effective masses are roughly equal to each other, i.e., the Fermi surface near half-filling 
is nearly circular. 

Consider now the electron-doped materials. Under the electron doping, the chemical 
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potential moves into the conduction band. The energy of a conduction fermion is E c = A + 
J(cos k x + cos ky) 2 — At' cos k x cos k y , i.e. it effectively has the sign of if reversed compared to 
the hole-doped materials. From the consideration above, we immediately conclude that the 
minimum of the electron dispersion is at (0, 7r) and the symmetry related points. The same is, 
indeed, true also for the selection due to quantum fluctuations in the nearest-neighbor model. 
Expansion around the minima yields two equivalent effective masses m± = ran = l/4|t'|. 

Finally, consider the electron dispersion in the Y — based hole-doped materials. The spin 
dynamics of the overdoped 123 systems was studied in a number of papers by K. Levin, Q. 
Si and coauthors 0. They found that to fit the photoemission data for Y Ba^Cu^O'j ||, one 
needs t ~ 0.25 — 0.3et> and relatively large next-nearest-neighbor hopping term t' ~ — 0.5£ 0. 
The values of the hopping integrals right near half-filling are not necessary the same as in 
Y BaiCuzO*; as the parameters of the effective one-band Hubbard model derived from the 
underlying three-band model generally depend on doping f27 |. We however simply assume 



that the values of t and t' change little with decreasing oxygen content. In this situation, 
\t'\ is very close to J. This implies that the mean- field hole dispersion, Eq(|(J), is nearly 
degenerate along k x = k y : the band minima at t = 0.3et> is at (vr,7r), but the quasiparticle 
energy at (vr/2, tt/2) is only 0.08et> above. To lift the near degeneracy, we calculated the 
leading self-energy correction to the hole dispersion in the expansion over the inverse number 
of orbitals. The procedure is described in some length in our earlier publication ||17|| , and we 
do not discuss it here. We found, that fluctuations stabilize the minima at (±7r/2, ±7r/2) up 
to much larger t' than in the mean-field theory. Specifically, the self-energy terms produce 
the energy difference, AE ~ 0.88 A (for 2S = 1), where AE = E^,*) — E^/i^/2)- For 
t = 0.3et> , we have A = At 2 < S z > / J ~ 0.9eu [^] and hence AE ~ 0.8et> . This implies 
that the actual critical value of \t'\ above which pockets are located at (ir, rr), is about 0.32et> , 
which is substantially larger than \t'\ ~ 0.13 — 0.15et> predicted for YBaCuO. 

The hole pockets at (±7r/2, ±7r/2) near half-filling are consistent with the results of 



photoemission studies of the insulating YBa 2 Cu 3 O e ^ @- These studies have detected 



some spectral features which can be interpreted as the dispersion through Fermi surface, 



but only close to the zone diagonal, i.e., near (tt/2,tt/2). Notice however that the closed 
Fermi surface near that point has not been restored experimentally. 

Another important point is that the critical value of t' is indeed model dependent - in the 
Hubbard-model calculations we found that it is larger than in YBaCuO, but, as we already 
mentioned in the Introduction, small cluster calculations for the t — t' — J model found the 
minimum of hole dispersion at (tt, 71") for the same values of parameters as we used ||. 

We further show how the different location of the pockets leads to a contrasting magnetic 
behavior near half-filling . 

C. Magnetic susceptibility 

In the SDW theory, the spin susceptibility is given by a ladder series of bubble diagrams 
(Fig.^). One fermion in the bubble should be above the Fermi surface, and one below. At 
half-filling, the only allowed combination is one fermion from the conduction and one the 
from valence band. Away from half-filling, the Fermi level moves into the valence band, and 
there are also bubbles with two valence fermions. The SDW expression for the susceptibility 
has been derived earlier [p!2j|l6| , so we quote only the result. In the static case, the total 
transverse susceptibility x + ~(<?) is given by 

1 - Uxo(q) 
where 

A 21 



Xo (l) = 2N ^ 



e k e k+q 



E d >\^\ L E k E k+q 



+ 



1 E 



Ef >| M | 



^ 6 k 6 k+q ^ 



Ek Ek+q 



,E k + ti k+q h k +h k+q/ 
1 



rpd rpd 

^k+q ~ ^k 



(10) 



Above we assumed that the system has a commensurate magnetic order. This requires 
that the static spin susceptibility be non-negative for all momenta (or, in other words, that 
all bosonic frequencies be real). Of special interest is the region near q = Q = (it, it) as 
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X 1 (q) turns to zero at (it, it) in accordance with the Goldstone theorem. Near this point, 
the static susceptibility has the form pOf 



X + ~ (q) = — ( 77^ 11 

Ps{q-QoY 

where N is sublattice magnetization (= 1/2 in our mean-field approach), and p s is the spin 

stiffness which should be positive. 

Let us first consider half-filling. Here only the first term contributes to xo- Performing 

an expansion in (|10D and substituting the result into (|9|), we obtain the "classical" spin wave 

result 



1 / 2(^ 
4 J £ 2 



^ = 7^(1-^- (12) 



Clearly then, the commensurate (it, it) state is stable at half-filling as long as \/2|t'| < t. 
This condition, though it may be modified by quantum fluctuations, is apparently satisfied 
in the LaSrCuO, NdCeCuO and YBaCuO families. 

We further consider the situation away from half- filling (5^0). Now we also have a 
contribution from the second term which involves only valence fermions. Expanding in this 
term around (it, it) and combining the result with flT2"|), we obtain 

p s (6) = p s (0) (1 - z) (13) 

where p s (5 = 0) is given by (0), and z is 

z = 4C/— hm > —5 —7 (14) 



>M 



k + q 



k+q k 



At small concentration of holes, the condition Ef < \p\ implies that the fermion with 
momentum k is within the hole pocket. For La and Y — based materials, these pockets are 
at (±7r/2, ±7r/2) where the sin 2 k factor in the numerator in ( |I4"D is approximately one. 
Accordingly, the summation over k yields the uniform Pauli susceptibility of free fermions, 
which in two spatial dimensions does not depend on the carrier concentration. Namely, for 
z we obtain 
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Z = 2C4 „S ( 15 ) 

TV 

For the case of (7r/2, tt/2) pockets, \t'\ < J, so that y^mlmji scales as 1/ J. In the mean-field 
theory, we also have C/ e // = £/ in which case z ~ U/J is a large number, and the spin 
stiffness immediately changes sign upon doping which means that the commensurate (tc, it) 
antiferromagnetic state becomes unstable. In more sophisticated calculations however, U e ff 
appears different from U because of the strong self-energy and vertex corrections in the large- 
U limit. In fact, the self-consistent solution for U e ff yields U e ff ~ J at U ^> t, and therefore 
z ~ 0(1) P,^l|,16|. Moreover, extremely close to half- filling, z vanishes logarithmically as 



z ~ 1/| log 5 1 because of the logarithmical singularity in the 2D scattering amplitude ||T6| . 
In any event, however, z rapidly grows with doping, and it is very likely that it quickly 
becomes larger than one in which case the commensurate antiferromagnetic state is no 
longer stable. Note that in the mean-field approach we are using, this instability does not 
imply a disordering transition, but rather a transformation into an incommensurate spin 
configuration. The equilibrium configuration at z > 1 has been discussed in our separate 
publication JT8[ . 



We now turn to the electron-doped systems. Here the hole pockets are formed around 
(0, 7r). Eq.(^) and fll4]) are still valid, but the numerator in (|14]) now vanishes right at the 
center of the pocket. Elementary calculations then show that because of the sin 2 k factor in 
(|l4l), z scales linearlywith doping concentration, and hence at small doping, p s acquires only 
a small correction 0(5). Clearly then, antiferromagnetism at (n, n) survives in the presence 
of a small density of electrons. This explains why N di- x C e x CuO ^ remains commensurate 
all the way down to the paramagnetic phase. 

We now discuss the width of the magnetically ordered phase. Within the present mean- 
field (or large S) approach, the on-site magnetization is nearly equal to its nominal value, 
and the rapid decrease in the stiffness in the hole doped 214 materials is not accompanied 
by the rapid decrease in the order parameter. In other words, the mean-field theory predicts 
that the system first becomes incommensurate and only then loses long-range order. There 

10 



are, however, numerous experimental reasons to believe that the disordering transitions in 
weakly doped high-T c materials are in the universality class of the nonlinear sigma model 
with the dynamical exponent z — 1 |32| (the most direct evidence is the observed linear 
behavior of the uniform susceptibility). This implies than in a more adequate model, the 
decrease in p s must eventually lead to the decrease in the sublattice magnetization such 



that both quantities vanish simultaneously. This is what has been found by Sachdev |33 
in the self-consistent large-N study of the Shraiman-Siggia model || in some range of the 
coupling constant values. In another range, he found an incommensurate transition within 
the ordered phase like in our approach. It is essential however, that the two scenarios differ 
primarily in the behavior of the sublattice magnetization with doping, while the doping 
dependence of the spin stiffness is nearly the same in both cases. In particular, for all values 
of the coupling constant in the Shraiman-Siggia model (where pockets are at (vr/2, vr/2)), the 
stiffness undergoes a rapid, nearly step-like, downturn renormalization under hole doping. 
We can, therefore, expect that the larger are the corrections to the stiffness at low doping 
(even if they are obtained in large-S* expansion, as in our approach) the smaller is the actual 
region of the magnetically ordered phase. Thus the width of the magnetically ordered 
phase in NdCeCuO should be much larger than in LaSrCuO. This is consistent with 
the experimental observation that magnetic order in NdCeCuO survives up to much larger 
doping concentrations than in LaSrCuO . A similar, though somewhat different, explanation 
of the difference of the magnetic phase diagrams of the two 214 compounds, based on the 
idea of localized electrons in NdCeCuO and mobile holes in LaSrCuO , was presented in ||. 
Finally, we notice that if the actual ratio of t'/t is such that the pockets in YBaCuO 
at low doping are located at (vr,7r), then the corrections to the stiffness scale linearly with 
x by exactly the same reasons as in the electron-doped compounds, and the commensurate 
(7r, 7r) configuration survives the hole doping. 
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D. Phase separation 

In SDW approach, we can also consider stability of the (it, it) phase at small doping 
against the formation of domain walls P4]JIo]j . This stability requires the longitudinal spin 
susceptibility to be positive. For the commensurate spin ordering, longitudinal spin fluctu- 
ations are always decoupled from transverse spin fluctuations, but at finite doping, they are 
coupled to charge fluctuations. The total static uniform susceptibility x zz can be obtained 
by straightforward manipulations starting from Eq.(52)-(54) in ||16| : 



o v Pauli 



y™ ~ 1± (16) 



where ^ Paub = ^mj_m\\/2ir is the Pauli-like susceptibility of doped carriers. If this sus- 
ceptibility is larger than 1/8J, the total longitudinal susceptibility becomes negative which 
signals the formation of domain walls. 

We calculated effective masses with self-energy corrections for all three types of materials, 
and found that for \t'\ ~ J/2, the tendency towards phase separation is nearly the same in 
the two 214 compounds (the Pauli susceptibility is slightly larger in LaSrCuO), and is 
much weaker in YBaCuO where the Pauli susceptibility is about two times smaller. At the 
same time, we found that the denominator in \ zz is positive in 214 materials, i.e., there is 
a stability against domain wall formation immediately away from half-filling. These results 
are consistent with the numerical analysis in Ref. ||. The latter paper also points to a 
possibility of a two-dimensional phase separation in the electron-doped materials which we 
didn't study. 

III. CONCLUSIONS 

To summarize, in this paper we presented simple calculations which show that the sta- 
bility of the commensurate antiferromagnetic state in the Mott-Hubbard insulators depends 
strongly on the location of the hole/electron pockets in the Brillouin zone. For LaSrCuO, we 
found pockets at (±7r/2, ±7r/2). The corrections to the antiferromagnetic spin stiffness from 
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the occupied hole states within these pockets rapidly grow with the carrier concentration 
and are likely to make stiffness negative, i.e., destroy commensurate spin ordering, already 
at a very small doping. On the other hand, in NdCeCuO, we found that mobile electrons 
form pockets at (0, tt) and the symmetry related points, in which case the corrections to 
the stiffness scale linearly with the density of carriers and do not destroy the commensurate 
spin ordering. We argued that the different behavior of stiffnesses is responsible for the 
experimentally observed difference in the widths of the magnetically ordered phases in the 
two 214 compounds. These results compliment the arguments and numerical analysis in Ref 

i- 

We also discussed the hole dispersion in YBaCuO and found that for the value of t' 
used to fit the photoemission data, band minima are likely to remain at (±7r/2, ±7r/2) 
though the quasiparticle energy at (tt, tt) is only slightly larger. This implies that the 
magnetic properties of LaSrCuO and YBaCuO are identical right near half-filling. At 
the first glance, this result seems strange as the hole pockets at (tt/2, it/2) apparently lead 
to incommensurability which has been observed in neutron scattering experiments only in 
LaSrCuO [§,0- However, these experiments were made only deep in the metallic phase 
when the Fermi surface is large and, in YBaCuO, is centered around (vr,7r). Recently, we 
considered |18| the situation when spin stiffness becomes negative while holes still occupy 
pockets at (7r/2, vr/2). We found that the equilibrium static spin configuration is not the 
planar spiral Shraiman-Siggia phase in which susceptibility is peaked at incommensurate 
momentum, but rather a noncoplanar configuration which very much resembles the (tt, tt) 
state and differs from it only in the existence of small transverse spiral component of the 
order parameter, S± ~ 0(y/x). For this configuration, the susceptibility still has a dominant 
peak at (7T, 7r). In other words, the spin structure adjustes to the negative stiffness of the 
(tt, tt) state in such a way that the peak position of the susceptibility does not change as 
long as hole occupy pockets at (tt/2,tt/2). Notice that this is consistent with the RPA-like 
analysis by Si et al in the metallic phase as very close to the magnetic transition they found 
the maximum in the susceptibility at (tt, tt) for both types of hole-doped materials. The 
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transformation of the hole Fermi surface with increasing doping content from small to large 
one, and the related change in magnetic susceptibility still need to be studied. 
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FIGURES 
FIG. 1. The RPA series for the total static transverse susceptibility. The first term represents 

the simple bubble, which is the building block of the ladder. Solid and dashed lines denote valence 

and conduction fermions, respectively. At half filling, only bubbles which contain one valence and 

one conduction fermion contribute to transverse susceptibility. 



